Quantum phase transition in spin glasses with multi-spin interactions 



On 



C3 



I 

c 

o 
o 



> 

O 

03 



c 

o 
o 



X 
S3 



Theo M. Nieuwenhuizen and Felix Ritort 

Van der Waals-Zeeeman Laboratorium and Institute for Theoretical Physics 
University of Amsterdam 
Valckenierstraat 65 
1018 XE Amsterdam (The Netherlands). 
E-Mail: nieuwenh@phys. uva.nl, ritortQphys. uva. nl 
(February 1, 2008) 

We examine the phase diagram of the p- interaction spin glass model in a transverse field. We 
consider a spherical version of the model and compare with results obtained in the Ising case. The 
analysis of the spherical model, with and without quantization, reveals a phase diagram very similar 
to that obtained in the Ising case. In particular, using the static approximation, reentrance is 
observed at low temperatures in both the quantum spherical and Ising models. This is an artifact 
of the approximation and disappears when the imaginary time dependence of the order parameter 
is taken into account. The resulting phase diagram is checked by accurate numerical investigation 
of the phase boundaries. 



I. INTRODUCTION 

The interplay between thermal and quantum effects in 
condensed matter physics is a longly debated problem 
The main differences between both type of effects 
relies on their dissipative nature. Thermal physics is in- 
herently dissipative and energy is not conserved while 
quantum physics is governed by the Schrodinger equa- 
tion where energy is conserved if the Hamiltonian does 
not depend on time. How to include relaxational effects 
in a systematic way in the regime where quantum effects 
are dominant is a very interesting open problem || . 

This question is of the most relevance concerning glassy 
systems (for instance structural glasses or spin glasses) 
which are manifestly non equilibrium phenomena. Re- 
cent developments in the understanding of the connec- 
tions between real glasses and spin glasses J§§ suggest 
that it is of interest to investigate that family of glassy 
models where the static phase transition is continuous 
from a thermodynamic point of view (i.e. there is no 
latent heat) but the order parameter is discontinuous at 
the transition temperature. These models are charac- 
terized by a one-step replica symmetry breaking (1RSB) 
solution at low temperatures and the existence of a 
dynamic singularity reminiscent of a spinodal instabil- 
ity Ji| . Let us summarize here the glassy scenario for 
this type of mean-field models. At a certain temperature 
(to be called Ta) in the paramagnetic regime the phase 
space splits up into different components or metastable 
states separated by high energy barriers (divergent with 
the size of the system), hence they have infinite lifetime 
in the thermodynamic limit. The number of these com- 
ponents is exponentially large with the size of the sys- 
tem A/" s = exp(S' c ) where S c is the configurational en- 
tropy or complexity. From a thermodynamic point of 
view the appearance of a large number of states does not 
induce a thermodynamic phase transition at Ta- Only 
at a "Kauzmann" temperature Tk lower than Ta a true 
thermodynamic phase transition (with replica symmetry 



breaking) is observed. At Tk the complexity S c vanishes. 
Hence, the glass transition is driven by a collapse of the 
complexity (entropy crisis) This is the mean-field 

version of the Gibbs-DiMarzio scenario for the glass 
transition. The dynamical behavior of the system in the 
region Tk < T < Ta is then dominated by the existence 
of a large number of components which trap the system 
for exponentially long time scales (t ~ e aN where N is 
the system size). Whether a sharp Tk exists in finite 
dimensions is still a largely unsolved problem (for recent 
numerical simulations sec JlT[| ) . 

In this direction, it has been recently proposed a ther- 
modynamic picture of cooling experiments in spherical 
p-spin models ]l2[ |. This new thermodynamic approach 
gives an explanation for the paradox of the Ehrcnfcst 
relations at the glass transition. The main new point in 
this approach is that the configurational entropy changes 
along the transition line. 

If the glass transition is driven by a collapse of the con- 
figurational entropy it is natural to ask how this scenario 
is modified in the presence of quantum fluctuations. Gen- 
erally speaking, quantum phase transitions appear when 
an external perturbation reaches a critical value at zero 
temperature. Because at zero temperature the entropy 
vanishes at any value of the external field it is expected 
that also the complexity should vanish everywhere at zero 
temperature (at least if there is no ground state degen- 
eracy, and this is the situation for the mean-field models 
we will consider here). In the absence of complexity it 
is natural to suppose that any adiabatic process at zero 
temperature (for instance, a process in which the exter- 
nal field is slowly turned off) could take the system to the 
ground state of the system. If complexity were not fully 
removed at zero temperature such expectation would fail 
since at zero temperature quantum tunneling processes 
could not take the system out of the traps during any adi- 
abatic process, mainly because the height of the barriers 
is extremely large p6fl . 

A hint to this problem was recently reported in H] 
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where it was shown that in a certain class of mean-field 
models where the Gibbs-DiMarzio scenario is valid, like 
the random orthogonal model ||, the complexity van- 
ishes at zero temperature. The transition turned out to 
be second order at zero temperature. That proof was ob- 
tained in the framework of the static approximation in- 
troduced by Bray and Moore How much general is 
this result beyond the static approximation and in other 
family of models (for instance, the quantum Potts model 
Q) is still unclear. 

A family of models which has received considerable at- 
tention during the past years are the spherical and 
Ising [^3) p-interactions spin-glass models. These models 
are also characterized by a classical continuous thermo- 
dynamic transition with a discontinuous jump in the or- 
der parameter. The purpose of this work is the study of 
the quantum phase transition in this family of models in 
an external transverse field. The Ising case has been al- 
ready considered in the literature and has revealed novel 
properties in the phase diagram. In particular, Gold- 
schmidt f25f| computed the phase diagram in the p — ► 00 
case, i.e. the random energy model of Derrida (hereafter 
referred as REM ^4j) in a transverse field. In this case 
the static approximation is exact and computations can 
be easily carried out. Goldschmidt found a phase di- 
agram with three different thermodynamic phases, two 
of them are paramagnetic and separated by a first or- 
der thermodynamic phase transition with latent heat. 
The existence of first order phase transitions in spin-glass 
models with a discontinuous transition has to be traced 
back to Mottishaw who studied the random energy model 
{REM) in an external anisotropy field |2(|. Computa- 
tions for finite p were done later on by Thirumalai and 
Dobrosavljevic p7| . They found that the thermodynamic 
first order transition line ended in a critical point. Such 
a critical point is pushed up to infinite temperature in 
the p — > 00 limit. Thirumalai and Dobrosavljevic went 
further and computed corrections to the static approxi- 
mation finding similar qualitative results at high temper- 
atures. Such an investigation has been recently extended 
by De Cesare et al. [^8| to the low T region. Corrections 
to the p — > 00 limit are generally complicated specially in 
the (3 —* 00 limit where the two limits have to be taken 
in the appropriate way. In a similar context, recent re- 
sults by Franz and Parisi j2i| also show the existence of a 
first order line when two replicas are coupled in the T — e 
plane where e is the strength of the coupling between the 
replica's. 

Some computations done in disordered quantum phase 
transitions involve the static approximation (hereafter re- 
ferred as SA) introduced by Bray and Moore [|l5| . This is 
a reasonable approximation close to the classical transi- 
tion line (in particular it predicts a decrease of the tran- 
sition temperature as the external field is switched on) 
but turns out to be inaccurate at low temperatures where 
dynamical correlations in imaginary time start to play a 



role (this is the reason why the approximation is called 
static). To clarify better the physical meaning of the 
SA we present an alternative derivation of the mean-field 
equations by introducing a solvable spherical version of 
the quantum Ising model. We will show that the same 
scenario for the quantum transition is valid in both the 
quantum spherical and quantum Ising models in the SA 
as well as beyond it. 

For pedagogical reasons we will analyze in detail first 
the spherical version of the classical model which is much 
simpler to solve. After that we consider a quantum ver- 
sion of the spherical spin system (recently introduced 
in (3l]]). We give all the details that occur in the def- 
inition and evaluation of the coherent state path inte- 
gral. After considering a few toy examples we analyze 
the quantum spherical p-spin model. 

Next we will consider the Ising case and estimate 
the corrections to the SA numerically solving the self- 
consistent mean-field equations. We will show that the 
approximation gives a reasonable estimate (within 10%) 
of the position of the line boundaries which progressively 
improves as p increases. Deficiencies of the SA for both 
models will be also identified at low enough tempera- 
tures, in particular reentrance of the T — V boundary 
line is observed. 

The paper is organized as follows. In section II we 
introduce the p-spin spherical spin glass model and a 
derivation of the thermodynamic behavior in a transverse 
field is obtained with classical and quantum spins. Sec- 
tion III presents the solution of the Ising case, the anal- 
ysis in the SA and also beyond it. Section V presents 
the conclusions. Finally some appendices are devoted to 
several technical points. 

II. SPHERICAL SPINS 

In this section we will consider multi-spin interaction 
spherical models without and with quantization. The 
spherical model is defined by 

^ = ~~ Jiiii—ipSi! Sf 2 ' ' ' &ip ~ r J]] Sf (1) 

ii<*2<---<ip i 

where T is the the transverse field. The indices 
*l)^2; • "ip run from 1 to JV where N is the number of 
sites. The Ji 1 i 2 ...i„ are couplings Gaussian distributed 
with zero mean and variance p\J 2 /(2N p ~ 1 ). The spins 
have m components and are subject to the spherical con- 
dition 

N m 

]T]T^ 2 = iVma (2) 

i— 1 a—1 

where a is a given constant. 

In what follows we first consider the classical case. 
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A. Classical situation 

The replica calculation for the classical spherical model 
without transverse field was described by Crisanti and 
Sommers (CS) |T(| . The steps are straightforward: 1) 
consider Z n ; 2) average it over disorder; 3) rewrite it 
terms of q a p = (I/AT) J2i=i Stasia'' 4) insert factors 



1 = / dq af} 5{q a p - — S L S ip) 

J — OO ^ „'_1 



OO plOO 

dq af 



Airi 



(3) 



-f- 1 = 



(qd - - £q) 



P 2 



-r 



iln^ 



= 



(8) 



(9) 



The latter equation expresses that dF/dx — 0. This 
means that we consider thermodynamic equilibrium. For 
a discussion of the thermodynamics of slow cooling ex- 
periments, see |Q. 

B. The paramagnet and its pre-freezing line 



A similar representation of the spherical constraints in- 
troduces as Lagrange multipliers the "chemical poten- 
tials" fi a . After these steps, one interchanges the order 
of integrals. The remaining integrals over Sf a are all 
Gaussian (this is the benefit of the spherical approxima- 
tion) and can be integrated out. One is left with an 
integral over q a p, q a p and /i, which can be taken at its 
saddle point. As for V — the fields q a a can be inte- 
grated out ]l6[ ] , and one ends up with the replicated free 
energy 



2/3F n = -log(Z») 



^g^-trlng 



(4) 



a/3 

/3r 2 

ma) h (m — 1) ln(/3/i a< 

Ma 



As for the case Y = we assume a one-step replica sym- 
metry breaking pattern. This involves parameters fi, qd, 
q, and x. These are the chemical potential (jj, a = /i), the 
replica self-overlap (q aa — qd), the overlap between dif- 
ferent replicas inside diagonal 1RSB blocks (q a p = q for 
(a, [3) inside a block, while vanishing outside the x x x 
blocks) and the breaking para-mater in the Parisi scheme 
(x is size of block), respectively. Note that qd are less 
than ma since the spins can turn perpendicular to the 
z-axis. Following CS we obtain for n — > the "classical" 
free energy F cl = F n /n 



/3 2 J 2 1 
.s.„c = y~ {<& ~ £,q P ) - - ln(qd - £q) - 1 



£ /3r 2 

- ln(q d - q) + f3fi{qd - ma) \-(m-l) ln(/?/i) (5) 

X [1 

where £ = 1 — x. Optimization with respect to p, qd, q, 
and x yields the saddle point relations 



T 2 (m-l)T 
qd H — n H = toct 



P (3J 2 



{qT 1 -q p - 1 ) + 



T 



(6) 



(7) 



Let us first consider the paramagnet, where q = 0. Like 
in the case of the p-spin Ising glass in a transverse field 
(see the next section), we find a first order transition 
line separating two paramagnetic phases. This line is 
comparable with the boiling line of a liquid and has a 
critical endpoint. To find it we insert (Q) with q = in 
eq. (^J) and obtain 



P/3J 2 p _i T 2 
(mer- q d )(— — + — ) 

^ qd 

-(m-l^qP-i + Tl) 

* qd 



(10) 



At large T and T this has just one real positive solu- 
tion < qd < ma. However, below a critical value 
T cep there is a regime of T-values where there occur 
three solutions rather than one. The outer ones are sta- 
ble, while the middle one is unstable. This critical end- 
point (cep) has coordinates (T cep , Y cep ), determined by 
dY/dqd — d 2 Y/dq 2 , — 0. From this point a first order 
transition line originates towards the spin glass phase and 
intersects it at the multi-critical point (T mC p,Y mcp ) p7| . 
Along this line there is a finite latent heat, that vanishes 
at the critical end point. It separates a small transverse 
field phase with large ordering in the z-direction (q^) 
from a phase with smaller (q^ ) ordering in the z-direction 
on the large field side. 

In analogy with wetting phenomena, where a pre- 
wctting line occurs off coexistence, we call this the pre- 
freezing line. In order to motivate this term, let us first 
explain the situation of first order wetting of a bulk fluid 
A by a thin layer of a fluid B |l^,|li|. At bulk coexis- 
tence of A and B phases there is a wetting temperature 
T w . For T below T w a finite layer ( "wetting layer" ) of B 
atoms will cover the A phase; for first order wetting this 
layer remains finite in the limit T — > T~ . For T > T w 
there will be an infinite B layer ("complete wetting"). 
When the fluids A and B are off coexistence there is a 
difference in chemical potential A/i. Let us take the con- 
vention that on the A/j, < 0-side the B layer is always 
finite. Then when A/x — > CP for T < T w , the B layer 
will reach its finite thickness discussed for Ap — 0. For 
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A/i > + it will be infinite however, leading to a dis- 
continuous transition. For T > T w , however, the thick- 
ness of the B-layer will diverge continuously in the limit 
A/i — > Q — , in order to be infinite at A/i = 0, and it will 
remain infinite for A/i > 0. In that temperature regime 
there is a continuous transition at A/i = 0. Thermo- 
dynamics requires coexistence at a first order transition 
line (called the "pre-wetting line") which separates the 
regime of continuous and first order wetting. The end- 
point of this line is called the "pre-wetting critical point" . 
The pre-wetting line, occurrence of hysteresis, and, near 
the pre-wetting critical point, scaling of the jump in cov- 
erage across the line have been observed in 4 Hc on Ce 
po| and for methanol-cyclohexane mixtures [|TJ 0] . 

In our spin glass a similar situation occurs. The SG- 
PM > line is a continuous transition line (in the sense 
that there is no latent heat), whereas we will find a finite 
latent heat at the SG-PM < transition. Also in this situa- 
tion a first order line with non- vanishing latent heat must 
emerge from the point (T mcp , T mcp ) and divide the para- 
magnet into two regimes. It is the line discussed, and 
by analogy we propose to call it the pre-freezing line. 
Its critical endpoint can then be called the "pre-freezing 
critical point" . 



C. The spin glass phase 

This discontinuity of the paramagnet has no analog in 
the spin glass. There is only one spin glass phase, namely 
the continuation of the paramagnet, with continuous 
qd at the transition line x = 1. We stress that this also 
holds also when the PM < phase is the thermodynami- 
cally stable phase: also then the (metastable) SG phase 
merges with the (metastable) PM > phase at the x = 1 
line. 

To check this continuity of the SG-phase, let us insert 
eq. (||) into eq. (||) and replace the x-dependence by 
dependence on a new variable rj via 



p - 1 - 77 q d -q 



Eq. 



then becomes 



In 



p-i (p- i-n)(r) + i) 



V 



prj 



(11) 



(12) 



which has a solution < r\ < 1. This shows that rj 
is independent of F and T. (For T = this was noted 
already by Crisanti and Sommers |l6| ) . Once 77 is known, 
we can choose x and solve q = (p — 1 — T))qd/(p — 1 — £77) 
from (|p). Eq. (g) will then yield 



qd = 



P-1-€V {2T 2 (p-l~r)f 



p — 1 — 77 \ x 2r qp{p — 1) 



(13) 



At fixed x we can vary T. We thus know qd and q, and 
therefore find a curve T(T). At small enough T two val- 
ues of x can lead to a given point (T, T) ; we need the 
smallest of these two rr-values. By varying x between 1 
and this procedure then uniquely determines the spin 
glass phase. 

The pre-freezing line intersects the PM > -SG transition 
line at a multi-critical point (T mcp ,T mcp ) Jl7| . 

Just as at r = the transition PM > -SG occurs with 
x = 1. This is a thermodynamic continuous phase tran- 
sition. The SG free energy exceeds the PM > one by an 
amount of order £ 2 . The transition PM < -SG is ther- 
modynamically first order and occurs with x < 1. The 
transition line is fixed by equating the free energies of the 
PM< and the SG solutions. 



D. Low temperature behavior 

The transition line between PM < and SG will continue 
down to T = 0. Everywhere along this line there will be a 
latent heat accorded by a jump in entropy. For studying 
the low T behavior we set 



q d = (1 + Tr)q 



which implies 

p- 1 -77 



rT: 



n 



p(p- 

2 V 



q p r 2 = 1 



It 



T) {l + rTJP- 1 -! 1 



qr (p — V)rT qr 
At T — this yields /i = (77 + l)/qr and then 
27?r 2 



9 - 



p(p-l)(7/+l) 2 J 2 
For p = 3 it can be solved exactly 



q 2 p = van 



7/r 2 



3(7/ + 1) 2 J 2 



(14) 



(15) 



(16) 



(17) 



(18) 



showing that the SG phase cannot exist at large T. For 
small r the spin glass phase is stable; for larger values it 
becomes meta-stable and for still larger values it will be 
unstable. The free energy may be expanded in powers of 
T 



F = F Q + T(-—\nT + F 1 ) + 0(T 2 ) (19) 



One finds 



F 



rj + 2 T 2 qr 
2pr 7/ + 1 



^1 = -(*?+ 1 + ( m - 1) ln(»7 + 1) + mlnqr) (20) 
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In the paramagnet PM < one has qd 
r /y/rrur, implying 



F = -i\ 

Fi = -(1 
2 V 



:) 



(21) 
(22) 



Equating the T — results we find a transition at some 
r c . For small T we get 



-Fsg - F PM = A(T - r c ) + ST 



(23) 



with A > because the stable phase has lowest free en- 
ergy, and 

1 



B 



-(?7-ln(?7 + l)) + 



m ma 
— In 

4 ma 



(24) 



2 4 toct — q 

also positive. For small T the transition line has a linear 
slope 



r, 



A' 



-T 



(25) 



showing that there occurs no reentrance. For the case 
m = 2 and in units where ma = 1 , the full phase diagram 
for p = 4 is shown in figure 1. In figure 2 we show 
the latent heat as a function of the temperature. It has 
been computed along the boundary lines starting from 
the critical point, following the PM < -PM > and the PM<- 
SG lines. Note the existence of a sharp maximum at 
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FIG. 2. Latent heat for the classical spherical model for 
p = 4 versus T along the boundary line which separates the 
PM < phase from the PM > (right part above the multi-critical 
point) and SG phases (left part below the multi-critical point). 
There is a maximum at the multi-critical point (indicated by 
the dashed line). 
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FIG. 1. Phase diagram for the classical spherical 
model for p — 4, ma = 1, m = 2. The 
multi-critical point and the critical point are given by 
T mcp = 0.3703, T mcp = 0.8208, T cep = 0.4767, Y cep = 0.5878 
while at T = 0, T c = 1.503 



E. Quantum spherical spins 

Due to the form ( |l9| ) the entropy of spherical mod- 
els diverges as (m/2)LuT for small T. The related zero 
temperature specific heat C = m/2, occurring due to 
the Gaussian nature of the spins, is analogous to the 
Dulong-Petit law of classical harmonic oscillators. In or- 
der to have a physical description in the low T regime, 
one of us recently proposed quantization by analogy with 
harmonic oscillators |3l). Here we present some details 
of this approach. It follows the standard Trotter-Suzuki 
approach of thermal field theories, see e.g. the book of 
Negele and Orland |§. 

The approach starts from the Trotter formula of the 
path integral representation of the partition sum, in 
which the coherent state representation of the identity 
is restricted to coherent states described by parame- 
ters which satisfy the spherical constraint. For a set 
of harmonic oscillators S = {Sf} with (i — 1, •••,JV, 
a = 1, • • • , to) a coherent state is defined by 

i S )= e s -^io>=n{f: i ^Ka>) (26) 

where Sl^ p is the creation operator of the harmonic os- 
cillator (i, a) and S" is a c-number. The coherent states 
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are overcomplete and have innerproduct 



(S'lS) =e s '' s 



(27) 



which can be checked in various ways. The coherent state 
representation of the identity in Fock space reads in gen- 
eral 



i,a K.nia—0 



-q dSf * dSf ^ s s 



e-» *\S)(S\ 



(28) 



It was proposed by one of us |3T| to enforce the spherical 
constraint by restricting this representation to coherent 
states which satisfy the spherical constraint 



N m 



S * ' S = S i* S i = Nm<7 



i—l a—1 

This is done by replacing 



^-spherical 



(29) 



(30) 



dSfdS? s .. Sl 



S)(S|5(S* - S-Nma) 



where the constant C will be fixed later. In a Trotter 
approach one calculates 



tr e 



-pH(S T 



tr 



M 



(31) 



with e = p/M. Between all factors exp(— eH) one inserts 
the coherent state representation of the identity. For our 
spherical spins we have to choose the truncated versions 
(p0[). Let us number them by j — 1,---,M. It was 
shown by Negele and Orland that for normal ordered 
Hamiltonians 

(S^le-^ 8 ^ 8 "')^-!) - e s;-s J _ 1 -eff(s;,s j _ 1 ) (32) 

The term exp(S* • Sj_x) arises from the overlap of the 
coherent states while the eH correction can be 

found by expanding the exponential, using SoplSj-i) = 
Sj_x|Sj-_x) and its Hermitean conjugate (Sj |S^ p = 
(Sj|S*, and re-exponentiating the result. Corrections are 
of order e 2 and can be neglected in the limit M — > oo (for 
a discussion, see [p2[). Introducing the imaginary time 
variable r = je = j/3/M and writing out the spherical 
constraints in terms of an imaginary valued chemical po- 
tential //(r), this leads to the coherent state path integral 
representation or thermal field theory for spherical spins 



Z = J DfiDS*DSexp(-A) 
with integration measure 



(33) 



[dS*DS = T[ r r d3(Sf(r))dSft(Sf(T)) 



involving the constant Cm — C to be fixed below and 
the action 

A = dr{S*(r) • + /j,(t)(S*(t) • S(r) - Nma) 



dr 

•ff(S*(r),S(r-dT))} 



(35) 



where dr = e and dS(r)/dr = (S(r) — S(r — dr)/dr in- 
volves S(r) due to eq. ( pp| ) and S(t — dr) due to eq. (|3^). 
The trace structure leads to periodic boundary conditions 
S(/3) = S(0), related to the bosonic nature of the spher- 
ical spins. One might be tempted to take the continuum 
limit of eq. (|33j) . However, note that this is a danger- 
ous limit since problems can arise that do not occur in 
our discrete formulation |3^] . For a concrete example we 
have discussed in Appendix A what is the origin of the 
problem. 



1. Free spins in a field 

The simple case of free spherical spins in an exter- 
nal field is already non-trivial pi] . This is because the 
spherical constraint couples the spins. Let us consider 
the Hamiltonian 



(36) 



We can introduce imaginary time Fourier transforms 
Si(r) =J2 S ^ e ~ lUT 

OJ 

S^ = ^$>(r) e — 



(37) 



where u> = 2-ktiT is a Matsubara frequency with 1 < n < 
M and r = jf3/M with 1 < j < M. Integrating out the 
spins we obtain 

Z = J Dficxp(-NA) (38) 

with (denoting imaginary times again by j = r/e), 

A = -ema ^ A*0'e) + m trj In B - e 2 T 2 ^ B jy ( 39 ) 

3 33' 



where 



B 33' = (! + ^U e )) S jj' - 6 3.3 



+ i 



(40) 



G 



with <5i.m+i = 1 due to the periodic boundary condition. 
We write 



(41) 



where \i is the saddle point value, that will turn out to be 
real, whereas the deviations fij are imaginary and turn 
out to be 0{N^ 1 / 2 ). We expand to second order in fij. 
The matrix 



B, 



1 



ep 



-B jy {n) 



(42) 



has diagonal elements 1 and off-diagonal elements —a 
— 1/(1 + ep). Its inverse is [||| 



i 



i-3 



,M A 



i > j 



1 



,M 



M+i-j 



I < J 



(43) 



We can now expand the action to second order in [ij. 
This gives after some algebra 



A = A + A 1 + A 2 
A is the saddle point free energy 

[3F = Aq = —/3fjtma + mln[(l + e/i) 



(44) 



M 



BY 2 
1 - — 
A* 



-Puma + mlnfe^ — 1] 



(3T 2 
ft 



This gives the saddle point equation 
m 



1 



-0H 



2 



At zero field one has = a/(a—l), yielding BF 
with infinite temperature entropy 



On 



m[a In a — [a — 1) ln(er — 1)] 



(45) 
(46) 

(47) 



Due to the scaling of the spherical constraint with m and 
the harmonic nature of the spherical spins, this yields an 
equal amount of entropy for each spin direction. 

For non-zero field the large temperature behavior is 
still of this form. For small temperatures, however, exci- 
tations will have a gap AE = [i(T = 0)=T / y / m(a — 1). 
This follows since, 



F 



(1 



-0AE 



(48) 



Note that this gap scales linearly in the field T, as 
expected for free spins in a field. Other quantization 
schemes have been proposed where the action involved 



a second oder derivative in imaginary time Q, [ p4| , 
p5| . Physically this is due to a kinetic term of the form 
(dS/dr) 2 rather than our first order derivative S*dS/dr 
arising from the Trotter approach. The kinetic terms de- 
scribe different physics, e.g. the quantized kinetic energy 
of a rotor. Such a system always has a finite energy gap 
due to its harmonic oscillator character. Spin systems 
are fundamentally different. Spins have no kinetic en- 
ergy, and for quantized spherical spins the energy gap 
indeed vanishes when the field vanishes. 

The next terms in eq. ( f44| ) fix the prefactor of the path 
integral. They are discussed in appendix B. 



2. Pair couplings 

Another non-trivial situation is quantized spherical 
spins that are coupled in pairs in the presence of an ex- 
ternal field. This covers both the ferromagnet and the 
spin glass cases. 



H— jjjSioySjop (s% ov + &j 



(49) 



We can diagonalize the coupling matrix, and introduce 
imaginary time Fourier transforms 



A _ — ILJT 



(50) 



where e\ is the normalized eigenvector of Jjj with eigen- 
value J\. Integrating out the spins we obtain the inten- 
sive free energy 



f3F = —/3fjLma + mM ln(l + e/i) 



(51) 



f3/j,ma + m / p(J\)dJ\\n(e 



- / p(J\)dJ\ 



(3(1 _ e /3J;0 



p- J\ 



(52) 



where T\ = £\ e i^i is the projection of the field on 
eigenstate A. The thermally averaged occupation num- 
bers are 



p - iQ u - J x e l 



(53) 



These results have been analyzed for a ferromagnet on a 
simple cubic lattice. One has J\ — > J(k) = 2J(cosk x + 
cosky + cosfc 2 ), with integration measure d 3 fc/(27r) 3 . 
In particular one finds a low temperature specific heat 
C ~ mT 3 ' 2 due to spin waves in m directions. Note 
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that the spherical constraint also allows longitudinal spin 
waves 

When spins are only coupled in the ^-direction, while 
the field acts in the transverse (x) direction, one has in 
the limit M — » oo 



@F = -0pma + 



d 3 k 
(2^)3 



+ (m - 1) ln(e^ - 1) - 



pr 2 



fi — 6 J 



(54) 



From these equations the zero-temperature quantum 
phase transition in a transverse field can be analyzed. 
Some results were given in ref . . 

The case of a spherical spin glass with pair couplings 
between the z-components will be useful for a check of 
the results of next section. Here Jy are random Gaussian 
with average zero and variance J 2 /N. The distribution 
of eigenvalues is the semi-spherical law 



P(J\) = 



1 



(55) 



with — 2 J < J\ < 2 J. For any M one can now calculate 
the thermal occupation numbers 

T 



p(J\)d,J x 



[i - iflu - J x e 1 ' 
2T 



H — ift u + ^{n-iVL^y -4J 2 e 2 * 



(56) 



Below the phase transition the system will have con- 
densed partly in the mode with largest eigenvalue 2J. 
Its occupation number Nq = S^jS^j will be extensive. 
For lo = one then has the expectation value 



q = qdu=o = q + 



The free energy reads 



2T 



+ - 4J 2 



(57) 



PF = -Puma + J p(J\)dJ x 

x ]Tln(l + e M -(l + e J A )e 2 ^ ) 

+ (m - 1) y ln(l + efi- e 2ltu ) - — 
+ p(fi - 2J)q 



(58) 



where the (m — l)-terms arise from the transverse spin 
components and the last term from the ordering field, 
respectively. Using 



/ p(Jx)dJ x ln(a - bJ x ) 



In 



a + \/a 2 - 4b 2 a - \/a 2 - 4b 2 



2{a + Va 2 - 46 2 ) 



(59) 



for a = 1 + e/i — e 2teuJ , b = ee 2zeuJ , we can express the 
integral as 

- 1 + Pp(q d -q)- HM(q d - q)) - ^-(q d - q) 2 + (60) 



^ [ -1 + P(p - i^)qdu ~ \n{Mq dLU ) 



P 2 J 2 2n 2 

— ^— e 



Variation of eq. wrt S 2J yields q = S^jS 2 j/N = 
or fi — 2J, q > 0. In the latter case eq. ( p7| ) yields 
= q + T/ J. We may therefore make the replacement 



2/3 Jg + /3 Jq d q 



This finally leads to 



P 2 J 2 2 



/3 2 J 2 



qd-q 



(61) 



^ = -^{fd - 9 2 ) - - 1 + M% - ma) 

^ qd — q 

BY 2 

- HMQd -q)-- + mln[(l + e M ) M - 1] 

A* 

+ + ^ ~ ^u)qdu - ln(/9(/i - if^)^) 



/3 2 J 2 



e 9cU 



(62) 



In next section we shall recover this expression as the 
p = 2, x -> case of eqs. @, @, @. The con- 
dition x = occurs due to absence of replica symme- 
try breaking. This should be expected since the system 
condenses in only one mode, the one having the largest 
eigenvalue p(| . When also short range ferromagnetic in- 
teractions are present, thermodynamics and correlation 
functions can be solved exactly. The largest mode may 
then be due to the onset of spin glass ordering or of fer- 
romagnetism. |f37f 

For low T the representation (|l|) shows that the spe- 
cific heat behaves as T 3 / 2 . One can also determine the 
time-dependent correlation function 



Qd 



qduje 



= I p(^A)dJ A — 



r( M -A) 



(-0 < T < 0) (63) 



At T — it is unity in r = and zero in + , while it 
decays as qd(j) ~ l 1 "! -3 ^ 2 for T —* —oo- 



F. Spin glass in a transverse field 

It is known that a given classical Hamiltonian may 
come from several quantum Hamiltonians. A similar 
situation occurs here. The simplest case is where the 
classical Hamiltonian contains complex valued spins p0[ , 
which can be replaced by operators. We thus consider 
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the case where H depends at each site either on the cre- 
ation or the annihilation operator. For p = 4 we have 
the Hamiltonian 



H 



-1 

IT 



N 

E 



J. ., , O z t Q z t CZ CZ 

u t]kl^>i op u j op u k op a l op 



i op I 



(64) 



with Hermitean Jijki |30[ |. (For odd p we have to add 
Hermitean conjugate terms; here they are included al- 
ready.) This means that for each quartet (i,j,k,l) with 
i < j, k < I there are four independent random vari- 
ables, J[ 2 and J{' 2 , each having average zero and vari- 
ance 9J 2 /N 3 . In terms of J1.2 = J[ 2 + iJ'12 the cou- 
plings in eq. (|||) read J ijk i = J% H j = J x + ij 2 and 
Jukj = Jkjii — Ji — *^2- These results hold for i < j, 
k < I. The other sectors follow, of course, by symmetry: 
Jijki = Jjiki = Jijik = Jjilk- A similar approach will 
work for general even p. 

In the replicated free energy we look for a saddle point 
with /x a (r) = n and fe Q (T,r') - (S*(t)*S*(t')> = 
q d (r - r') and with <Z Q/3 (r,r') = (^(t)*S|(t')) = q a (i 
independent of r, t' for a ^ (3. For M Trotter steps 
the free energy for a one step replica symmetry breaking 
solution with plateau q and breakpoint x = 1 — £ reads 



F — Fstatic + F t 



quant 



(65) 



with 



B 2 J 2 ,„ p p , 1 



2 x qd - q 

BT 2 

+ B\i{q d - ma) - ]n[8fi(q d - q)\ - — + m ln[e^ - 1] 

/' 

being mainly twice as large as F c i ass i C in (H) , due to dou- 
bling of spin degrees of freedom (now complex, previously 
real). A more important difference is the replacement 
raln/3/z — > mh^e^ — 1). As we shall see, this improves 
quite a bit on the not-too-low temperature behavior. Af- 
ter deriving this expression at finite M we have replaced 
a term mln[(l + e/i) M — 1] (see also eq. (|62|)) by its limit 



,ln[, 



1]; we shall come back to this point below. 



The quantum correction reads at finite M 

P Fq Uan i 

= / J-l + S(n - iCl u )q dul - ln(8(p - ifl^qdu)] 

P 2 J 2 P/ 2 f 1 ^ P/2/ \ , P 2 J 2 ~p ian\ 

~ 1^ q d ( £ + T > q d ( £ - T > + ~^~ q d ( 6? ) 

r 

Here we have Fourier transforms 



qd 



qd^e qdc 



= h^ qd{T)e ^ r (68) 



and denoted q~ d = <idu=o- 

If one inserts in the J 2 term of F quant that qdij) = q d 
is independent of t, the J 2 -terms cancel. Then the 
g w 's can be solved, after which the whole F quant van- 
ishes identically. This is closely related to the static 
approximation (SA) of Ising models (see below) intro- 
duced by Bray and Moore |l5| ]. This approximation ne- 
glects the time dependence of the correlator q d {r) where 
%(t) qd being also independent of r. The remain- 
ing difference with previous classical theory is the re- 
placement mln/3/z — > mln(exp(/?^) — 1) in going from 
2F c i ass i C to F s tatic- This replacement already improves 
the low temperature behavior. The phase diagram in the 
SA can be numerically computed and is shown in fig- 
ure 3. It is qualitatively similar to that computed in the 
classical case (see figure 1). We find a thermodynami- 
cally first-order transition line with a multi-critical point 
terminating in the high-T phase in a critical end point. 
At low temperatures the first order line shows the phe- 
nomena of reentrance and negative latent heat along that 
line down to T — 0. This is a failure of the SA as we will 
show below. 

Beyond the static approximation the saddle point 
equations read: 



4-° 



F 2 ^ 

— 7T + (m — 1 —5 

\i 2 K ' - 1 

n — <f 



Ad ( e - T )<id ( £ + T ) 



4- =0 



2M 
pp 2 J 2 



l"- 1 = T-~ 



= 0^ - 



B 2 J 2 



{qd - q){qd - £q) 

£,q q 



x{q d - C?) 



(69) 
(70) 

(71) 
(72) 



The latter equation is solved by x = (p — 1 — r)) idd — 
q)/(rjq) with the same rj as in the classical case eq. (|12|) . 

Both for u =/= and for lo = we multiply 
dF quant /dqduj by g^w and go to the time domain. This 
yields 



8 Wd (r) + M (q d (T) - q d (T + e)) 



(73) 



MS Tj3 - — 



{qd - q){qd - £q) 



P (3 2 J 2 
2M 



J^qdir + e-T'yJ^ie-T'yfie + T') 



with t = je, j = 1, 2, • • • , M, and e = /3/Af . Note that 
now eq. ( |70|) becomes redundant, as it follows already by 
summing (|73[) over r. 
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1. Numerical solution at finite M 

We have numerically studied the quantum equations 
with m = a = 2. In order to compare with the classical 
case where (1/A) ^ Sf — 1, we have rescalcd spins S — > 
S/s/ma, S* — > S* /y/ma to yield a unit constraint also in 
the quantum case. This amounts to scaling temperatures 
and fields to T — > T/(mcr) 2 , T -> r/(mcr) 3 / 2 . 

The Trotter limit M — > oo should be taken. This set 
of non-linear equations can be numerically solved for a 
different values of M. We find that depending on the 
regularisation term [ j38| we use for the static part of the 
free energy eq. (|6^) the low temperature behavior of the 
first order line shows quite strong finite M corrections 
and the numerical extrapolation to the limit M — > oo is 
not safe. To overcome this problem we did the follow- 
ing: we took two different regularisations for F sta tic, he. 
we replaced the term mlog(e' 3M — 1) by the general mi 
dependent expression, 

(m - m x ) log(e^ - 1) + m x log((l + ^-) M - 1) (74) 

Note that in the limit M — > oo this expression coin- 
cides with the term mlog(e^ — 1) for any value of mi. 
For finite M the behavior of the first order transition line 
at very low temperatures strongly depends on the value 
of M. This is a direct consequence of the correct order 
of limits in the saddle point equations where the limit 
M — > oo should be taken before the limit T — ► 0. Con- 
sequently the behavior of the line in the limit T — > is 
quite different if these limits are taken in the opposite 
way (i.e first, T — > and later on M — > oo). 

In the classical model we had M = 1, mi = m. We 
found that the first-order line matches the T = axis 
with a negative slope without reentrance (see figure 1). 
In the quantum case with mi = 1, M > 1 this situation 
persist with smaller value for r at T = and larger slope 
of the transition line. In the static approximation (eq. 
(|66f)) one has mi = 0, so the regularisation term coincides 
with the M — > oo limit itself. In this case the model 
shows the phenomena of reentrance for finite M close to 
zero temperature (like in the SA of the Ising model, sec 
below). As we expect the transition line to have infinite 
slope, this indicates that for each M an optimal value for 
mi exists between and 1 where the slope is infinite. 

Numerically we proceed in the following way: for a 
given value of M we determined the value of mi such that 
the first order transition line meets the T = axis with 
infinite slope. The optimal lines obtained in this way for 
different values of M (we took M = 1, 3, 5, 7, 9, 11, 21, 41) 
are then extrapolated to M — * oo. We found that a sec- 
ond degree polynomial in 1 /M is enough for such extrap- 
olation to be accurate (even if higher order polynomials 
are needed at very low temperatures). 

At zero temperature we obtain T c — 0.5453 for p = 4 
and mi ~ 0.033. The resulting extrapolated boundary 



line for p — 4 is depicted in figure 3. Only at very low 
temperatures deviations form the SA are important. The 
phenomena of reentrance has now dissapeared since this 
was an artifact of the SA. The latent heat at very low 
temperatures, across the first order transition line, van- 
ishes exponentially with 1/T. At high temperatures cor- 
rections to the SA are indeed very small and the value 
of the transition field in the SA is always larger than the 
exact M — * oo extrapolated value. The opposite result 
is found at very low temperatures. 

In section III we will see that a similar scenario is valid 
for Ising spins in the SA and also beyond it. 
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FIG. 3. Phase diagram for the quantum spherical 
model for p = 4. The multicritical point and the crit- 
ical point are given in the static approximation (SA) by 



0.486, Tc, 



0.2118, r c 



0.3743. 



The continuous line is the SA and the dotted line is the 
M — > oo extrapolation which yields F c = 0.5453 at T = 0. 



2. Continuum limit: M — » oo 

Let us now take the limit M — ► oo. We set 

qd(T)=q + p(T) (75) 

with p(r) = p(r + j3). Eq. ( |73| ) implies a discontinuity 
for r = + : p(0 + ) = p(0) — 1 with p being left-continuous 
at 0. At other r one gets 



dp(r) 



^-j^ +p ^C dT ' HT) - piT+T ' )] 

x [(q+pir^-^q+pi-r'^-qP- 1 } (76) 
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Further one has 



r 2 

q + p(o) + — 



p(p-l)J 2 



v p-i 



(/drp(r))^ 



(77) 



(78) 



(79) 



2 /-/3/2 



2 - nP- 1 ] 



P 4- I ' dr'[(g + Kr , )) p/2 - 1 (g + K-^)) p/ --'/ 

L J-/3/2 

The internal energy reads 

[/ = - / dr[( g +p(r)f/ 2 (g + p(-r))f/ 2 -^] 

J-0/2 

.7-/3/2 M 



(80) 



For p = 4 its Fourier representation reads 

U = —f3J 2 2_j QdunQcUj 2 Qdui3Q(lui4 

(p-i-^gP-i 2 r 2 

Pf> w =o 

77 



(81) 



where = Plu + c/<L,o- 

These equations are particularly useful at T = 0. It 
can then be seen that p(r) ~ 1/t 2 for r — > ±00, implying 
that Pui = T J&T P (T)e iuJT ~ T(l + \u\ + iuj). Expanding 
the sums in ( pl[ ) in powers of p we find sums over 1,2, and 
3 frequencies. The one-frequency sum can be calculated 
as follows. We extend the Euler-Maclauren formula to 
complex functions with non-analyticities of the form |w|, 
and obtain 



„ 1 du> 27rT 2 m d/, 



j=2t™T 



6 dui 



(82) 



The two and three frequency sums are convolutions of 
this and produce also T° + T 2 terms. This yields a be- 
havior U = Uq + U2T 2 , which implies a linear specific 
heat, C w 2U2T, in the spin glass phase. This result 
holds for the present p-spin model. Unlike stated previ- 
ously Q, C ~ T also holds for the p = 2 + 4 model 
with infinite replica symmetry breaking. For p = 2 only 
there is no replica symmetry breaking, and the system 
is in another universality class. As discussed above, one 
then has C ~ T 3 / 2 . 



III. ISING SPINS 

The Ising p spin-glass model in a transverse-field is 
defined by, 



n 



il<i2<---<i P 



(83) 



where erf , erf 



are the Pauli spin matrices and V is the 
the transverse field. The indices ii^is,- "ip run from 
1 to N where N is the number of sites. The Ji x i 2 ..-i p 
are couplings Gaussian distributed with zero mean and 
variance p\J 2 / (2N P ~ 1 ). We shall choose units in which 
J = 1. 

Here we will compute in the SA the phase diagram of 
the model eq. (p3|) and show that coincides in its essentials 
with that reported in the previous sections. Detailed 
computations of the quantum Ising model eq.(p3[) have 
been already presented in the literature. Here we only 
sketch the main steps of the derivation of the saddle point 
equations skipping the details. The interested reader will 
find more details about their derivation in f25[ , ^7| , fj8| . 

The free energy of the model is computed using the 
replica method as in the previous section. After discretiz- 
ing the imaginary time direction using the Trotter-Suzuki 
decomposition we obtain a problem described by an ef- 
fective Hamiltonian, 

H eff =AJ2 Jy °\°] + C ( 84 ) 

i<j t it 

where the time index t runs from 1 to M and the 
spins a\ take the values ±1. The constants A, B and 
C are given by A = jj;B = ± ln(coth(^)); C = 
ln(| sinh( 2 ^-)). Now we apply the replica trick and 
compute the average over the disorder of the replicated 
partition function, 



/n 
[dJ] £ exp^ Ht ff 



(85) 



Computations are easily done and the problem can be 
reduced to a dynamical equation involving Ising spins in 
a one dimensional chain. The free energy reads, 



rwO n 



(86) 



where, 



mA) -¥ + i Tr(QA) 



A 



abtt' 



(87) 



with Q" fc , A" fc being the order parameter and the trace 
Tr is done over the replica and time indices. The term 
if (A) is given by, 

H(A) = £ exp(]T ^ £ K^alai +3^ 

a ab tt' at 
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The most general time traslation invariant solution of 
these equations is given by, 

Qai=Q ab (a^b); A«'=A ah (a + b) (89) 



Qaa = qd(t~t') 



A*=\ d (t-t') (90) 



Because at zero transverse field the classical solution 
is a one step of replica symmetry breaking we look also 
for solutions of this type in the quantum case. We divide 
the n replicas into n/m boxes K of size m such that m 
divides n. The saddle point solution when a ^ b takes the 
form Q% = q; A"' = A if a, b e K and Q% = A«' = 
otherwise. If a = b we have Q** a = q d (t — t'),A** a = 
X d (t — t'). Finally the free energy reads, 



P 2 (m - 1) 



2 9X+^J2 qd (t-t')Kt-t') 



1 r°° 

— In / d Px E m (x)) (91) 
m l 



and dp x = dxexp(— x 2 )/(2tt)^ is the Gaussian measure. 
The order parameters are determined by solving the sad- 
dle point equations, 



dj_ = dj_ = di_ 

dq dX dm 
df df 



= 0; 
= 0; 



dq d {t-t>) d\ d (t-t>) 
where H(x) is given by, 

S(z) = ^exp(e(x,{a t })) 

{°t} 



(92) 
(93) 

(94) 



with 



A 



e(x,{a t })) = —J2(Ht-t')-X)a t a t , 



AVXx^ 



(95) 



Solving equations (B2h,rt93|) we get, 



A = |g p - 1 ; X d {t-t') = ^{ qd {t-t')Y- 1 :, (96) 

q=«{alf»] q d {t-t')=«aTa r F» (97) 
(l-p)(3 2 q p _ 1 
rn 2 



= — ln(/ dp x Z m (x))- 



1 



<< ln(S(x)) >> (98) 



where the averages << .. >> and (.) are defined by, 



where 0(x, {<Tt}) is given in equation (p5|). 

The solution of this system of coupled equations is 
quite complex because there is an infinity of parameters 
(q d (t—t')) which needs to be computed in a self-consistent 
way. For p = 2 (the SK model in a transverse field) the 
transition is continuous in the presence of the transverse 
field and there is only one quantum paramagnetic phase. 
For p — 2 these equations have been studied using five 
different methods. These are: 1) doing a self-consistent 
approach [|?J] or a Ginzburg-Landau expansion |4(J, 2) 
Performing exact small M calculations |yj 3) Perturba- 
tive expansions in the field J42] , 4) Numerically solv- 
ing the Schrodinger equation p43| and 5) doing quantum 
Monte Carlo calculations p4|j45| ]. In the case p > 3 the 
transition is discontinuous and eqs.(|97j),(98) have been 
perturbatively solved by expanding around the p — > oo 
limit p7p^ ] where the SA (see below) is exact. Here we 
will revisit the SA showing that the phase diagram of 
the model coincides in its essentials with that presented 
previously for the spherical model. We will go beyond 
the SA later on and numerically solve eqs.(|9^),(p7|),([9^) 
by doing finite M calculations in order to check the reli- 
ability of that approximation. 



A. Zeroth order solution: the static approximation 

The SA amounts to consider q d {t) and X d (t) indepen- 
dent of t. This corresponds to supress quantum fluc- 
tuations. This is exact at zero transverse field but it 
turns out to be inaccurate at finite field and crucial for 
the thermodynamic properties at zero temperature. The 
failure of the SA is very clear in case of continuous quan- 
tum phase transitions where the quantum critical point is 
characterized by the dynamical exponent z, an exponent 
which cannot be computed within the SA. The situation 
is slightly better in first order quantum phase transitions 
where there is no critical point. Hence there is no diver- 
gent correlation length and imaginary time correlation 
functions can be well approximated by constant values 
p6[ . Generally this approximation can be the source of 
pathologies at low temperatures where the third principle 
of thermodynamics is usually violated. At not too low 
temperatures we will see that this approximation yields a 
phase diagram in qualitative and quantitative agreement 
(within 10 per cent in the worst case p = 3) with the full 
dynamical solution. 

Let us first analyze the solution of the mean-field equa- 
tions in this approximation and study the phase diagram 
of the model. 
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Introducing A<j(t — t') = Xd, qa(t — t') = qd in the free 
energy eq.(|9l|) we get 



0/ 



4 



— In / dp B (S(a?)) ro (101) 
to 7 



where, 



/oo 
dp z cosh(T(x, z)) 
-oo 

T(x,z) = (6 2 + /3 2 r 2 )2 ; 



(102) 
(103) 

The A, Xd, m arc determined by solving previous equa- 
tions eqs.©,©,© and q, qd are determined by solv- 
ing the following equations " ~ ] 



(sinh(T)^) 2 » (104) 



,,=<<<co 8 h(T)<A, 2 + <2I!fffi)>> 



(105) 



where the average < < . . > > was previously denned in 
eq.(p^) and the average (..) is given by, 



B(x,z) 



E(x) 



(106) 



The phase diagram of the model can now be computed 
for any arbitrary value of p. As in the spherical case we 
find two different paramagnetic phases. Putting q = X = 
and to = in eqs.(104),(105) they reduce to a single 
equation, 



qd = 



J-oo ^ X $(x) 

JZo d P* cosh($(a;)) 



(107) 



with *(a?) = V/3 2 r 2 + P 2 ^dx 2 and X d = pq^ / 2 - 
This equation can be numerically solved. Like in the 
spherical case one finds two paramagnetic solutions sep- 
arated by a first order transition line with latent heat. 
Let us call QP > and QP K the quantum paramagnetic 
phases associated to the largest and smaller value of qd 
respectively. The transition line can be constructed us- 
ing the Maxwell rule. As temperature increases the la- 
tent heat decreases. Consequently the first order line 
ends in a critical point (j3 c ,T c ) with mean- field critical 
exponents. The existence of this critical point has been 
already pointed out by Dobrosavljevic and Thirumalai 
p7[ . We have numerically computed it for several val- 
ues of p. Details of these computations are given in the 
Appendix C. In the infinite p limit this critical point is 
pushed up to infinite temperature pq ] and its scaling be- 
havior in the large p limit has been analitically obtained 
in J2l finding T c = 0.7579, T c = 0.2593^. 



As temperature is lowered the first order line finishes 
in a multicritical point which separates the three phases 
of the model (two paramagnetic QP > and QP < and one 
quantum glass QG ). The boundary lines which separate 
the paramagnetic phases from the QG phase correspond 
to different thermodynamic phase transitions. The line 
which separates QP > from QG has no latent heat (this 
is the continuation of the usual first order classical phase 
transition at V — 0) . The line which separates QP K from 
QG is a first order transition with latent heat. The latent 
heat is positive when crossing the QP < — ► QP > line as 
well as the QP < — > QG line. Lowering the temperature 
the QP < — ► QG line is determined by the Maxwell con- 
struction but allowing to to be different from 1 and q, X 
jump to a higher value when crossing the QP < — > QG 
line. For low values of T the breaking point m is nearly 
proportional to the temperature and the difference be- 
tween qd and q proportional to the temperature in the 
paramagnetic (in this case q is equal to and qd is propor- 
tional to T) as well as in the quantum glass side (where 
q and qd reach a finite value smaller than 1). 

The behavior of the latent heat in the boundary lines 
QP < — > QP > and QP K — ► QG as a function as a func- 
tion of the temperature is the following: starting from the 
critical point (where there is no latent heat) and lowering 
the temperature the latent heat increases as a function of 
the temperature reaching a maximum in the multicritical 
point. Then the latent heat decreases and vanishes like 
Tip- 1 ) a t i ow temperatures. 

We have analyzed in detail the phase diagram for two 
different values of p. We have chosen a small (p = 3) 
and a large value of p (p = 10). The phase diagram for 
p = 3 is shown in figure 4 and that of p — 10 is depicted 
in figure 5. The latent heat along the thermodynamic 
first order transition line is shown in figure 6 (for p = 
3, 10,oo respectively). For sake of completeness we also 
show the dynamical transition line for different values of 
r in the QP phase (see jl4| for more details how this line 
has been computed in the random ortoghonal (ROM) 
model). The main result concerning this dynamical line 
is that it crosses the first order transition QP > — QP < 
below the ending critical point. 

The main difference between figures 4,5 is that for p = 
3 the critical point is hardly observable (but it is there!, 
the difference between T cp and T mcp being of order 10 -5 ). 
Also the latent heat corresponding to p — 3 is smaller 
than that of p = 10. Being the case p = 3 so close 
to p = 2 (where the transition is continuous and there 
is no multicritical point) it is natural to find that the 
transition is nearly continuous. Note that also for p = 3 
the dynamical and the static transition lines are both 
very close to the multicritical point. 
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line has associated a dynamical transition line (the dy- 
namical transition predicted in the framework of Mode 
Coupling theories) which meets the static line precisely 
at the multicritical point. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

r 



FIG. 4. Phase diagram for p = 3 with Ising spins. The 
critical point is given by T cp = 0.3528, T cp — 1.1078. The 
multicritical point is extremely close to the critical one and is 
indistinguishable from it in the figure. At zero temperature, 
F c = 1.174. The dashed line is the dynamical transition. 

As anticipated in the previous sections we observe in 
figure 6 that the latent heat becomes negative at very 
low temperatures. For p — 3 this happens below T ~ 0.1 
while for p = 10 this effect persits but is hardly observ- 
able. This is a small effect because the latent heat is 
already of order — 10 -3 for p — 3 and — 10~ 5 for p = 10. 
The same comments presented in the spherical model 
also apply here. A negative latent heat implies reen- 
trace close to zero temperature. Consider the Clapeyron 
equation for first order transition lines ^ — L/(TAM X ) 
where L is the latent heat and AM X is the change in 
transverse magnetisation when crossing the PG < — * QG 
line. Because AM X is always negative (increasing T the 
transverse ordering M x increases) a negative latent heat 
implies ^ > 0, i.e. reentrance. In fact, reentrance is 
observed in figures 3 for p — 3 and hardly observable 
(but there is) in figure 4 for p = 10. In the limit p — > oo 
reentrance dissappears |2j|. Like in the case of spherical 
quantums spins reentrance for finite p is an artifact of 
the SA. 

Perturbing around p — 2 we expect the following sce- 
nario to be valid: for p = 2 the transition is continu- 
ous (there is no latent heat) and there is no multicriti- 
cal point. Above a critical value p^p > 2 it appears a 
multicritical point which separates a first order transi- 
tion line (with latent heat) from a thermodynamic sec- 
ond order transition line. The second order transition 
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FIG. 5. Phase diagram for p = 10 with Ising spins. 
The critical point and multicritical point are given by 
T cp = 0.7765, T cp = 0.8903, T mcp = 0.5543, T mcp = 0.854. 
At zero temperature r c — 0.8855. The dashed line is the 
dynamical transition. 

In the regime 2 < p < p^ there is a unique quantum 

(2) 

paramagnetic phase. Above a given value p c such that 
Pc 2 "* > Pc a first order transition line appears with two 
paramagnetic phases in both sides. Wether p c is larger 
or smaller than 3 is unclear. Within the SA, we expect 
Pc to be quite close to 3. A definitive answer to this 
question requires a full analysis of the theory beyond the 
SA. In this sense a perturbative study in p = 2 + e would 

(2) 

be useful. The fact that p = 3 is close to p c explains 
why the transitions looks like a continuous one with very 
small latent heat (see figure 6). 

B. Beyond the static approximation 

As said in the previous section it is natural to expect 
that the SA works well enough if the transition is not 
continuous. In fact, we expect it should yield better and 
better results when p increases (in the p — > oo limit it 
is exact) even if it is always wrong because it definitely 
violates the third law of thermodynamics at zero tem- 
perature p7[ . For smaller values of p it should be pro- 
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gressively worse being uncontrolled close to the quantum 

... (2) 

transition point at p ~ p c . 



0.18 
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0.08 
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of the two phases we can obtain the transition point for 
different values of M. Then we extrapolate the free en- 
ergies, latent heat as well as the transition point, to the 
M — > oo limit. A second degree polynomial in 1/M fits 
quite well the data. 

In figures 7 and 8 we show the free energy as a func- 
tion of the transverse field T for p — 3,T = 0.3 and 
p = 10, T = 0.4. The M — ► oo extrapolation is com- 
pared to the static ansatz which appears to be a reason- 
able approximation in this case. The error in predict- 
ing the value of the critical field is ~ 10% for p = 3 
(yextrap = 1.021 ± 0.002, Tf 4 = 1.14604) and 2% for 
p = 10 (T e c xtra P = 0.841 ± 0.001, Tf 4 = 0.8798). This 
error should increase at lower temperatures. The latent 
heat is shown in figure 9 for different values of M as 
well as the extrapolation to M ~ * oo compared to the 
value obtained in the SA. The agreement is very good 
for p — 10 but not for p = 3 where the SA predicts a 
latent heat nearly 4 times larger than expected. 
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FIG. 6. Latent heat for the p-spin Ising spin glass model 
with p = 3, 10, oo (from left to ri ght) . The last case are the 
results obtained by Goldschmidt p5|. It is shown along the 
boundary lines which separate the PM K phase from the other 
phases as a function of the temperature. There is a maximum 
at the multicritical point. 

To go beyond the SA we have numerically solved 
eqs.©,©,© for different values of M for a fixed 
value of /? and extrapolating the results to the M — > oo 
limit. This is a method which usually yields good re- 
sults and has been applied in several cases to continuous 
quantum phase transtions in disordered systems Ej^E^ . 
The essentials of the method has been already presented 
in section III.F for the spherical quantum model. Here 
we will show how the method works for first order quan- 
tum phase transitions in Ising models. Our procedure 
is quite simple: we solve the system of non linear equa- 
tions (|96"l),(|97|),(|98l) for different values of M looking for 
a quantum paramagnetic QP < and a quantum glass QG 
solution. We have used periodic boundary conditions 
such that <7m+i — T The QP < solution is described 
by q = A = and q^ P {t — t') different from zero. With- 
out much effort the equations can be solved in the QP < 
phase up to M ~ 16. In the QG phase the solution of the 
set of non-linear equations requires more computational 
effort (because q, A and m are now finite and some one 
dimensional integrals cannot be avoided). In this case 
we were able to solve the equations only up to M = 12. 
Looking at the crossing point between the free energies 
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FIG. 7. Free energy as a function of the transverse field V 
in the Ising case for p = 3 and T = 0.3 for different values 
of M. The dotted lines correspond (from below to above) 
to M = 4, 6, 8, 12. The long-dashed line is the free energy 
extrapolated to the M — * oo limit. The dashed line which 
connects the filled circles contains the transition points for 
different values of M, the last one is the extrapolated tran- 
sition in the M — > oo limit. The continuous line is the free 
energy in the SA and the star indicates the transition point 
in that approximation. 



15 



0.12 



-0.9 - 



-1.0 - 



-1.1 - 



-1.3 - 



-1.4 



-1.5 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



0.10 



0.08 



a 

X 



1-J 



0.04 



0.02 



0.00 

0.00 0.05 



0.10 0.15 0.20 0.25 

1/M 



FIG. 8. Free energy as a function of the transverse field F 
in the Ising case for p — 10 and T — 0.4 for different values 
of M. The dotted lines correspond (from below to above) to 
M = 2,4, 6, 8, 10, 14. The long-dashed line is the free energy 
extrapolated to the M — > oo limit. The dashed line which 
connects the filled circles contains the transition points for 
different values of M, the last one is the extrapolated tran- 
sition in the M — > oo limit. The continuous line is the free 
energy in the SA and the star indicates the transition point 
in that approximation. 

Another interesting result in figures 7 and 8 concerns 
the jump in the transverse magnetisation. Using the re- 
lation Ai x = — §p- this jump manifests in a discontinuous 
change of the slope of the free energy as a function of T. 
From the figures it can be observed that the transverse 
magnetisation always decreases going from the QP < to 
the QG phase. The jump is very small for p = 3 and 
increases for larger values of p. 

It is very difficult to perform numerical calculations 
at much low temperatures, mainly because the scaling 
behavior in M is found when the ratio -p- is small in 

order to extrapolate to the continuum limit jjy — ► 0. 
At T = 0.1 we have studied the case p = 3 for values 
M = 8,9,10,12,13,14, this last case being the limit of 
our computational capabilities. The results are shown in 
figure 9 where we plot the latent heat as a function of jj . 
It is difficult to extrapolate to M — > oo because we do 
not have large enough values of M in order to do that. 
The data is compatible with the fact that at very low 
temperatures the latent heat is negligible in the M — > oo 
limit. 



FIG. 9. Latent heat for the Ising case, p — 3 (filled 
circles) and p = 10 (filled triangles) at the transition 
point for T — 0.3, T — 0.4 respectively as a function of 
1/M. The dashed lines are second degree polynomial fits 
in 1/M to the data. The triangle (p = 10) and the cir- 
cle (p = 3) in the vertical axis are the values estimated 



in the SA (L 



r extrap 
^p=3 



extrap 
p=10 



0.027 ± 0.001, L 

SA 



SA 
p=W 



0.02767, 



= 0.0097 ± 0.001, L£f 3 = 0.036) 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have investigated the quantum phase 
transition in spin glasses with multi-spin interactions in a 
transverse field. We have introduced a solvable spherical 
model which yields a phase diagram qualitatively similar 
to that found in the Ising case (in the static approxima- 
tion -SA- and also beyond). Details of the quantization of 
the spherical model have been given. We find indications 
that the specific heat is linear at low T. This is possibly 
related to a finite density of two level systems in the free 
energy landscape. We have also seen that p-spin models 
in a transverse field (spherical and Ising) typically have 
a first order transition line in the paramagnet, that we 
have called the pre-freezing line. 

For the p-spin models the study indicates (as expected) 
that the static approximation (SA) can be considered as a 
classical approximation where quantum fluctuations are 
fully neglected. A zero order calculation shows that the 
SA seems to yield qualitative good results for first order 
(but not too weak) phase transitions at not too low tem- 
peratures. The situation is different for continuous quan- 
tum phase transitions. In particular we have checked 
the approximate validity of the SA in both the spherical 
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(with quantized spins) and the Ising model numerically 
computing the free energy and the transition line. This 
has been done solving the time correlator q^(t — t') for 
finite values of M and extrapolating to M — > oo. The SA 
predicts the phase diagram of the model with reasonable 
accuracy. For instance, for p = 3 in the Ising case the SA 
yields the phase boundaries with a precision within 10% 
improving for larger values of p. The approximate valid- 
ity of the SA is restricted to high temperatures. Indeed, 
at very low temperatures the SA fails. This manifests in 
the phase diagram of both the Ising and spherical cases 
(this last one with quantized spins) which display the 
phenomena of reeentrance. This pathology is related to 
the incorrectness of the SA and disappears when taking 
into account quantum fluctuations. 



the dependence (if any) of the dynamical exponent z with 
e. Recent results in the ROM model suggest that 
the quantum dynamical exponent could be not univer- 
sal within mean-field theory. This suggests that models 
with the same classical behavior may display different 
quantum behavior in presence of the same type of per- 
turbation. 

It would be very interesting to investigate the problem 
of the existence of more than one quantum paramagnetic 
phase in the quantum Potts model where it has been 
suggested (like in the ROM model) that the transition 
becomes continuous at zero temperature p8[ . These are 
subjects for future research. 

Acknowledgments. F.R is grateful to the Founda- 
tion for Fundamental Research of Matter (FOM) in The 
Netherlands for financial support through contract num- 
ber FOM-67596. The authors acknowledge hospitality at 
the ISI (Turin, Italy), where part of the work was done. 



0.080 



0.060 - 



*j 0.040 - 



0.020 - 



0.000 1 1 1 1 1 ' 1 ' 1 

0.000 0.040 0.080 0.120 0.160 

1/M 

FIG. 10. Latent heat for the Ising case, p — 3 at the transi- 
tion point for T = 0.1 as a function of 1/M. Extrapolation to 
the M — > oo limit is not safe since we are far from the scaling 
region. 

In the simplest scenario the multicritical point should 
appear as soon as 1 < pc < P- In this case the phase 
diagram should be qualitatively similar to that of figure 
4 with only one quantum paramagnetic phase. Above 

(2) 

a second critical value p c the multicritical point would 
develop a line ending in a critical point restoring the two 
different quantum paramagnetic phases like is observed 
in figure 5. It would be interesting to understand (in the 
spherical as well as in the Ising cases) how the phase dia- 
gram of the model changes when expanding in p = 2 + e. 
If 2 < p < p^ the transition should remain continuous 
for small e. Then it would be interesting to investigate 



APPENDIX A: OTHER DISCRETIZATIONS OF 
THE COHERENT STATE PATH INTEGRAL 

The coherent state path integral has an obvious expres- 
sion in the continuum limit M — > oo, dr — > 0. However, 
that is a dangerous limit, which may introduce problems 
that do not occur in its finite M expression f32| . A typi- 
cal case is the following sum over Matsubara frequencies 
lu = 2%nT (n = 1,---,M) 

P = ^ln(/3/z- iflD-jy (Al) 

where A/" is an appropiate normalization and 

ft w = iMT(l - e iew ) w w (M < 1) (A2) 

This sum can be carried out after expanding in powers 
of e leuJ and yields in the limit M — ► oo 

P = ln[(l + en) M - 1] -> ln[e^ - 1] (A3) 

provided we choose JV = M In M. The common ap- 
proach, however, is to approximate Q u ps u>, choose M = 
and to extend to sum from —M/2 < n < M/2 
— ► — oo < n < oo, which yields the result 

P = In 2 sinh -0fi = ln[e /3AI - 1] - 0(i/2 (A4) 

This ill-defined procedure thus brings a different result 
for the non-singular part. Those terms also show up in 
the zero point energy, that is to say, terms that may arise 
when normal-ordering of the creation and annihilation 
operators. The common approach also yields a different 
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answer for the first derivative of P wrt /i. For the second 
derivative the convergence is quick enough to yield the 
same answer in both approaches. 



APPENDIX B: NORMALIZATION OF THE 
PATH INTEGRAL: FREE SPHERICAL SPINS IN 
A FIELD 



The second term in eq. (E3) is 



Ai = {—ma - 



1 - a 



AI 



r 2 



(Bl) 



As expected, it vanishes when /i is taken at the saddle 
point. The quadratic fluctuations yield 



„ . 9 \ ^ o , m r . 

2 ^ = ^ 2 E^(t^m + ^) 



(B2) 



2 \ - , ma M 

jj' 



r 2 alJ'-'"l +a M -l^"l 



1 -a 



,1/ 



) 



The /x 2 and fijfiji terms can be calculated by going to 
Fourier space. Using the equation of motion Ai can be 
rewritten as 



2A 2 = Me 2 ^ MwM-o>[<7 + <L,o 



1 



/j, 2 1 — ae 4 
Thus the ^-integrals yield 

Cm 



1 — ae 



(e^ - l) 2 
E7)] (B3) 



Z 



2M^ 

/3 M 



,-Ao-iD 



(B4) 



D = V In ; aK . w ; — B5 

^ (1 - ae lW£ )(l - ae" !we ) 

The w-sums can be carried out using 

ln(l - foe ±lwe ) = ln(l - b M ) (B6) 

For the leading behavior at T ^ we get with b = 1 — 
D » Mini +21n(l _ « VmJ^- 



If we choose 



C M = (27rMiVCT) M / 2 



(B7) 



(B8) 



it follows for r 7^ that the free energy has, on top of the 
extensive part TA$, a non-universal contribution of order 
7V°M 1 / 2 . For r stricktly equal to 0, there is a universal 
term N° In M. Both terms are non-extensive and can be 
omitted if one first takes N large and then M . Actu- 
ally this is also the limit that underlies the saddle point 
approximation. Physically it is also the natural limit, as 
for fixed small T the M = oo limit is reached already for 
M ~ 1/T independent of N. This example shows that 
the extensive part of the free energy of quantum spheri- 
cal spins is a well defined, natural object. Non-extensive 
parts are more delicate. 



APPENDIX C: CRITICAL ENDPOINT IN THE 
STATIC APPROXIMATION: ISING CASE 



Here we give the equations which yields the critical 
endpoint in the SA in the Ising case (see also [^7) for the 
original derivation). 



Starting from eq.(107) we define the function g(x), 

(CI) 

pqP- 1 /2 and 



9(9) = Q - ®(x)~ 1 x 2 sinh<i>(x) 



where $(x) = \/(3 2 T 2 + /3 2 Xx 2 and A 
the average (.) is defined by, 



A{x) 



SToo dp* cosh(x) 



(C2) 



The paramagnetic phases are found by solving the 
equation g(qd) = 0. This yields one solution at very high 
temperatures and three solutions at lower temperatures. 
Of these three solutions two of them are stable (the ones 
with largest and smallest values of q£) while the other 
one (that with an intermediate value of qj) is unstable. 
This is the same scenario as in the spherical model. The 
critical point is then determined by the coalescence of 
these two stable solutions. This gives the equations, 



9^ d ) = (^) q=qd = (^) q=qd =U 



These three equations read, 



(5 2 p{p - 1) 2 (2) 

4 Id fc , (Qd) = l 



(C3) 

(C4) 
(C5) 



(2-p)fP(q d ) = 



PMp-i^- 1 



/ C (3) M (C6) 



where /j , /j , /i are the first three cumulants as- 
sociated to the functions f( n ' (n = 1, 2, 3), 



/W = $-1^2 snln <j) 
/ (2) = $" 3 a; 4 ($cosh$ 

/( 3) 



sinh $) 



(C7) 
(C8) 



$~ 5 x e ($ 2 sinh $ - 3$ cosh $ + 3 sinh $) (C9) 
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These equations can be exactly solved yielding 
T cp , T cp , <f? for different values of p. Is not difficult to 
generalize this set of equations beyond the SA in the 
general case. 



APPENDIX D: EQUATIONS FOR THE ENERGY 
IN THE ISING CASE 



In this appendix we give the exact expressions for the 
internal energy used to compute the latent heat in sec- 
tion IV. B. We start from equation ([nj) by evaluating the 
derivative u = ^§q-- This yields, 



dC_ _ f3{m - l)qP _ J3_ 
d(3 ~ 2 



2AP 



£(?<«(* -op 



p(m - l)q\ + ^ £ ~ *W - f ') 



where C and S(x) were defined in eqs.(Q) and ( |94| ) 
respectively and dp x is the Gaussian measure. Doing the 
last integral by parts and rearranging terms we get the 
final expression, 



2/3T f3(m-l)qP 
M = _rcoth(— ) 



2M w smh(-^- 



(D2) 



In the continuum limit M — > oo the <7d(i) becomes a 
continuous function of time yielding, 



dqa 
dt Jt=o 



(3{m - l)qP 1 f p 



(q d (t)) p dt (D3) 



In the QP < phase at zero temperature in the large 
r regime we have qd{t) ~ exp(— tT) yielding u ~ — T — 
l/(2pr). Note that the SA is only exact in the limit 
p — > oo where the energy is given by it = — V . 

It is also easy to check that in the SA the energy is 
simply given by, 



?(«S-(i 



m)q H 



/?r 2 « — v ; yy >> 



T(x, z) 



(D4) 



where the averages (..), << (..) >> and T(x,z) where 
previously defined in eqs.(B9f), (lOq), (103) respectively. 
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